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1. Introduction 

Investigation of various aspects of gravitational waves is one of tlie main topics of 
contemporary relativistic physics. In rigorous treatments of gravitational radiation in 
the framework of Einstein's theory, most studies of this problem have concentrated 
on radiative spacetimes which are, at least in some directions, asymptotically flat. 
Until recently, exact "cosmological" gravitational waves have attracted less attention. 
One of the reasons for this is that the definition and characterization of gravitational 
radiation in such a more general setting is not obvious. 

Indeed, standard techniques used for investigation of asymptotically flat 
spacetimes (such as the Bondi-Sachs approach which allows one to introduce the 
Bondi mass and momentum, and to characterise the time evolution and energy balance 
including gravitational radiation in terms of the news functions — the analogue of the 
radiative part of the Poynting vector known from electrodynamics) cannot be applied 
when the cosmological constant A is non- vanishing. Moreover, even in the relatively 
simple case of vacuum spacetimes, but with A ^ 0, it had long been known [1-3] 
that the radiative component of the gravitational field is not unique. It depends 
substantially on the null direction along which it approaches a given point at conformal 
infinity I because it has a spacelike character when A > or timelike character when 
A < 0. 

It was recently demonstrated in the series of papers [4-11] that such a directional 
structure of gravitational or electromagnetic radiation can be generally described in 
a closed form. In fact, it has a universal character that is essentially determined by 
the algebraic type of the given field, namely by local degeneracy and orientation of 
its principal null directions on I. The general behaviour for both the cases A > 
and A < was, respectively, described in [8] and [9], and subsequently summarized 
and presented in a unified way in the review [10]. However, more specific and explicit 
analyses are highly desirable to exemplify these results. To that end, it is necessary 
to consider suitable classes of exact radiative solutions of Einstein's equations which 
admit a nontrivial cosmological constant. A particularly interesting family of such 
spacetimes is the so called C-metric which can be physically interpreted as describing 
the field of (possibly charged) non-rotating black holes that miiformly accelerate in 
the de Sitter or anti-de Sitter universe. Detailed studies of the directional structure 
of radiation and related properties of these spacetimes have already been presented 
in [6,7]. 

The purpose of the present work is to extend and further generalize these results, 
namely to consider an effect of a twist. In particular, we will consider the large 
Plebahski Demiahski family of exact type D solutions of the Einstein- Maxwell field 
equations [12] (see also [13]). This not only includes the famous Kerr-Newman rotating 
black holes, the Taub-NUT spacetime, and the (anti-)dc; Sitter metric, but also their 
arbitrary combination. In addition, the resulting black holes may accelerate due to 
conical singularities (such as cosmic strings) attached to them. The general form of 
the metric thus contains seven free parameters which characterize the mass, charge 
(electric and magnetic), rotation, the NUT parameter, acceleration of the sources, and 
the cosmological constant. 

This Plebanski-Demiahski class of solutions has recently been described in more 
detail in [14-17], where a list of further references can also be found. In particular, 
a new form of the Plebanski-Demianski metric has been presented in which the 
parameters were given a clear physical interpretation and from which the various 
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well-known special subfamilies can easily be obtained. Two new parameters directly 
representing the acceleration of the sources and the twist of the repeated principal 
null congruences have been introduced, and then (following the idea first applied in 
the particular C-metric cases [18,19]) the remaining freedom has been used to simplify 
the roots of the corresponding metric functions. This has helped significantly in the 
interpretation of this family of solutions: in particular to distinguish the Kerr-like 
rotation and the NUT-like contribution to the twist. A simple unified form of the 
previously known black hole spacetimes has thus also been obtained. We will employ 
its straightforward coordinate modification in the present analysis. 

This paper is organized as follows. In section 2 we will review the Plebanski- 
Demiahski family of solutions, introducing convenient coordinates and frames. In 
section 3 we will explicitly derive the specific radiative field components of the 
Plebahski-Demiahski spacetimes, using the approach given in [11] which applies to 
general type D metrics. We will cover all possible cases which may occur, both for a 
positive and negative cosmological constant. In the final section 4 we will analyze and 
visualize the thus obtained directional structure of radiation and, as a main result, its 
amplitude in various points at conformal infinity. 



2. The general Plebaiiski Demianski family of black holes 

In [12] Plebahski and Demianski presented a complete class of type D solutions of 
Einstein's equations, with a possibly non-zero cosmological constant A, which may be 
either vacuum or electrovacuum (such that both repeated principal null congruences 
of the Weyl tensor and that of the non-null electromagnetic field are aligned). As 
shown in [14], it is useful to write the general Plebahski-Demiahski metric as 



ds^ = 7^ ( - ujp'^daf + „ „ „ {ujdr + r'^daf 

(1 — apr)"^ \ -|- u'^p^ r^ + w^p^ 



-^d^^ + ^^d.^ 



(1) 



where V{p) and Q{r) are quartic functions 
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Q = (uj^k + + g^) - 2mr + er^ - 2a— r^ - (a^k , 

LU \ 3 



(2) 



and m, n, e, g, A, e, k,a and uj are arbitrary real parameters of which two can be 
chosen for conveniencelj Depending on the specific choice of these parameters (and 
thus the number and types of the roots of the polynomials V and Q) and on the 
admitted range of the coordinates, the above class of solutions contains a very large 
number of various spacetimes with different geometrical and physical properties. 

Physically most interesting are those solutions which describe black holes. It was 
demonstrated in [14-17] that, in such a case, it is particularly convenient to express 
the parameters e, n and fc, which occur in the metric functions ([5]), in terms of new 
parameters a and I as 

Lj^k I ,9 9, / , 9 9 9, A \ 

4a-m-{a^ + 3f) —{u;^k + e^ + g^) + - , (3) 



l'^ — P LO \ LO'^ '3 



X The original metric given in [12] is recovered by setting a = 1 = oj and putting k = -y — — A/6. 
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(5) 



1 + 2a-m - 3a^ — (e^ + q^) - PA 

UJ UJ^ 

The remaining scaling freedom can be used to set the twist parameter lu to any 
convenient value, provided a and / do not both vanish. For the particular case 
a = = I we apply the relations e — 1 — a^{e'^ + g^) , n — , k — 1. 

We are thus left with seven physical parameters, namely A, m, e,g,a, Z and a. 
These have direct interpretation as representing the cosmological constant, mass 
parameter, electric charge, magnetic charge, Kerr-like rotation, NUT parameter and 
the acceleration of the black holes, respectively. Indeed, performing the simple 
transformation 

I a „ , {I 



P - 
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the Plebahski-Demiahski metric ([T]), ^ becomes 



where 
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Q, = l (/ + acos^^r , 

UJ 

p2 = ,'2 + + acos6')2 , 

P ~ \ — a^ cos 9 ~ 04 cos^ 9 , 

Q = (w^fe + + 5^) - 2mr + e 



(8) 
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with e, n and A; in Q given by ([3])- ([5]). It is also assumed that jaaj and |a4| are 
sufRciently small that P has no roots within the considered range 9 G [0, tt]. 

This general metric immediately reduces to the Kerr-Newman-NUT-de Sitter 
solution when the acceleration vanishes, a — 0, and further to the familiar forms 
of the Kerr-Newman-de Sitter black hole spacetimes when the NUT parameter / is 
set to zero or to the charged NUT-de Sitter spacetime when the rotational Kerr-like 
parameter a is set to zero. Then, putting the charges e, g and the cosmological constant 
A to zero, the well-known forms of Kerr (and Schwarzschild) and Taub-NUT metrics 
are recovered. On the other hand, when a ^ but we set I = and then a = w = 0, 
the non-twisting C-metric solution is immediately obtained which describes charged 
black holes uniformly accelerating in Minkowski, de Sitter or anti-de Sitter universe 
(according to the sign of A). Therefore, the line element (O is a very convenient metric 
representation of the complete class of accelerating, rotating and charged black holes 
of the Plebahski-Demiahski class. 

Let us also recall that the spacetime has a Kerr-like ring singularity at p = when 
\a\ > \l\. When \a\ < \l\, the metric is singularity free. Although these two cases have 
identical metric forms ([7]), their singularity and global structures differ substantially. 
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2.1. Alternative form of the metric 

For an investigation of the asymptotic behaviour of gravitational and electromagnetic 
radiation it is necessary to cover completely the regions near the conformal infinity T. 
However, for the Plebahski-Demiahski family of space-times this is difficult to achieve 
using the coordinates of the above metric ^ . The main reason is that the coordinate 
r e (0, +oo) does not extend up to X , which is located at = , in all directions, that 
is for all values of 6. It will thus be convenient to introduce another metric form of the 
Plebahski-Demiahski family of black holes. This is obtained from ([T]) by introducing 
the reciprocal radial coordinate q = —r^^. The resulting metric[§| reads 



d.^ - ^ (- -| (dr - Vda)2 + ^ (u^q^dr + daf + ^ dp' + ^ dq' ] , (10) 
where 



n = -{q + ap), 
0^ = 1+ Lj'^p'q' . 



Vip) = k + 2—p - + 2amp^ - a'Ujj'^k + e' + g 



2/, ,2;, , „2 , „2\ , ^^^2 



A 

3' 



(11) 



A 



Q{q) = - ( a^fc + ^ ) + 2a-q + eq' + 2mq^ + {ijj'k + + g') q^ 

^ 3 / LU 



The parameters e, n and k are still given by ([S])-®. The roots of Q given by Q{q) = 
represent horizons at specific values of g < . To retain a Lorentzian signature of the 
metric (|10l) it is necessary that V > 0. The coordinate p must thus be restricted 
to a particular range between appropriate roots P = which correspond to poles of 
the black holes. This is automatically achieved by assuming the parameterization 
p = {l/uo) + {a/ijj) cosd with 9 e [0, tt]. In fact, to express all Plebahski-Demiahski 
black hole spacetimes — including those with a non-trivial NUT parameter I — in a 
closed and explicit form it is useful to finally apply the specific linear transformation 
([U on the metric (fTO| . This would yield the metric analogous to ([7]), only with the 
"radial" coordinate r replaced by q via the above relation q = — r^^. Notice also that, 
consequently, > in the whole spacetime. 

In order to characterize the directional structure and asymptotic behaviour of 
radiation, it is first necessary to express the Weyl and the electromagnetic tensor in the 
natural frame of null vectors k, 1, m, m, such that gap gap m°'m'^ = -1-1 . 

We need to distinguish two possible situations, in which the region near the conformal 
infinity T is either stationary (Q > 0) or non-stationary (Q < 0). 

In the regions when Q > 0, the null frame is 

" ^ ^ {dr-ioq'd,) + VQdg^ , 



V2g \VQ 
fl ( 1 



idr-Loq'd,)-VQd, 



rdp + ^iup'dr + d,)] , (12) 



V2gV ' VV 



Vdp-^iup'dr+d^) 



V2gV ' VV 



§ In fact, for a = I = lo this reduces to the original form of the Plebaiiski- Demianski metric given 
in [12], see (2.1) and (3.3) therein. 
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VVdp-^{u;p^dr + d„)] , (13) 



The only non-zero component of the Weyl tensor in both these frames is 



*9 



(m + in) + (e^ + g^) 



1 + i ujpq 



ap 



1 — i ujpq 



(14) 



which explicitly demonstrates that the spacetimes studied are of algebraic type D 
and that the (future oriented) tetrad vectors k and 1 are aligned with the double 
degenerate principal null directions of the Weyl tensor. The only non-zero component 
of the Ricci tensor associated with the aligned electromagnetic field is 



4 



Z Q 

which shows that for e ~ Q ~ g the spacetimes are vacuum. In fact, the corresponding 
vector potential is 

^ _ e g ( dr - ujp^da) -gp{da + tjgMr) 

~ T~, 2 2 2 ' 

1 + uj^p'^q^ 

so that the electromagnetic field F = dA is given by 

F Q-^^[e[l-uj'^p^q^) - 2gujpq] dg A dr - [g(l - ^V?^) + 2eujpq] Ap h Act (17) 

— ujp'^[e{\ ~ uj^p^q^) — 2gujpq\ dq A dcr — ULiq^[g{l — uj^p^q^) + 2eujpq\ Ap A dr| . 
With respect to the null frame (|12p we thus obtain 

$1 = ^|e(l-a;V?')-25WTO-i[g(l~c^V'?') + 2ec^TO]|, (18) 

while with respect to the frame (fT3|) we obtain the complex conjugate of (fT8|) . In both 
these cases the relation $ii = 2$i$i yields (fTS]) . 



3. Asymptotic structure of radiation 

We will now investigate the behaviour of radiation in the above family of solutions. 
3.1. The conformal infinity 

The first step is to identify the conformal infinity X of the spacetimes. It turns out 
that the above metric pU|) is particularly convenient because it is naturally related to 
a suitable conformal (unphysical) metric gap through the conformal factor > 0, 

gap = Vl'^gaf3 , Vl = -{q + ap). (19) 
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The conformal infinity J is located at $7 = which obviously occurs at 

I : g = -ap, (20) 

where the conformal metric gap evaluated on scri, that is 5^/5 |i, is regular. Since p is 
finite, X occurs at finite q, so that the conformal infinity of all Plebanski-Demiahski 
black hole space-times is completely covered. In particular, in the equatorial plane 
given by p — it is simply located at g 0. This parameterisation is also well behaved 
in the limit with vanishing acceleration (a = 0), in which case the whole scri is given 
by g = 0, that is r = 00. This is a standard result for Kerr-Newman-(anti-)de Sitter 
black holes. 

It is then necessary to find the vector n normal to conformal infinity X. This 
is defined as fi" = N g"^ dfjfl where TV is a suitable function, cf. [10]. Since the 
gradient of is df} = ~{dq + adp), it follows that n -Ng^^{Q dq + aV dp). To 
set the appropriate normalization gajsn^n^ = —sign A on conformal infinity, it is 
necessary to choose = -pj . Indeed, the vector 

n = -^^g-HQdq + aVdp) (21) 

normal to X has the norm gapn°'n^ = jj-^ Q~^{Q + cP'V) , in which the fimction 

Q + a'V = -^(1 + a'cjV) - ^ + + 9^){q " ap){q^ + aV) (22) 

+ 2m{q^ — aqp + a^p^) + e{q — ap) + 2a[n/uj) , 

evaluated on scri (where 17 = 0) obviously reads 

A , , .X A 

3' 

Therefore, on X it follows that [gai3n"n'^)^ = —1 or +1 when A > or A < 0, 
respectively. Thus, in the asymptotically de Sitter spacetimes n is the unit timelike 
vector, while in the asymptotically anti-de Sitter spacetimes it is spacelike, as required. 
The relation ^ also implies Qx = -f - {V{p) + |wV)j • Moreover, 

V{p) > Q everywhere. (24) 

For A > the function Q is thus negative at all points on I, which means that the 
complete region close to the conformal infinity is non-stationary. In the case A < 
the situation is more complicated: the regions where Q > are stationary while the 
alternative regions where Q < are non-stationary, and these are separated by the 
Killing horizons located at Q = 0. This structure occurs also at the conformal infinity 
(which has the anti-de Sitter- like character). For small values of the acceleration 
parameter a all regions of X are stationary while for large acceleration the non- 
stationary region are also present: 

Qi > iff (V{p) + |c.V) < -| , (25) 

Qi<0 iff a' (Vip) + |c.V) > -| • (26) 

In the equatorial section where p = 0, these conditions simply reduce to a^k < 
■|, respectively. 



[Q + a^V)^ ^-^{1 + a W) - -iA ■ (23) 
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3. 2. The structure of radiation when A > 

We will first derive the explicit form of the field components which describe the 
radiation in the Plebanski-Demianski family of spacetimes in the case when the 
cosmological constant is positive. In such a case the conformal infinity 2 has a 
de Sitter-like character with the unit normal vector ([2T|l being timclike. 

The function Q is everywhere negative near and on such scri, so that we will 
employ the null frame k, 1, m, m given by (fT3|) . which is adapted to both principal 



null directions. Since gaf3 k^n^ ~ ~y 7i \l ^ — 9aP l^n^^ , where n = n, the 
frame is future-oriented and it also satisfies the normalization condition ei — €2, see 
equations (18) and (19) in [11]. The associated orthonormal tetrad ts,qs,rs,Ss: 

^"72(*^ + ^^)' l=73(ts-qs), ^{r,-is,), (27) 

adapted to these algebraically special null directions, is thus 

n , — 



Q 

^ -1 



^ ^ (28) 



Q 

Ss = ^ [UJP Or+ da) . 

Now, it is natural to choose the following reference orthonormal tetrad adapted to X: 



:{dr- Ujq da) , 



(29) 



So = 1= ( Ujp dr+ Oa) ■ 

g Vr 

Using (|23p . (pij) it is obvious that to = n = iln on X, i.e. to is a future-oriented 
timelike normal to T. Also, qs = qo , Sg = Sq , and 

ts = cos^^^s to + tan6's To, rg = cos~"'^0s To + tan^s to, (30) 

where 



sin6's = y'"^^^' cos6's = yi - q;2^^. (31) 

(These are explicit relations of those given generally in [11], see equation (21) therein.) 
Consequently, the magnitude of the radiative components of the gravitational and 
electromagnetic fields are described by the expressions (27), (28) therein, namely 

l*4l = ^^^A{e, CP, 6^) , = L\3^A{9, CP, e.) , (32) 

|?7| cos'^b's r/^ cos^Ps 

where 

A{e,4>,es) = (sine + sin cos 0)^ + sin^ cos^e* sin^ . (33) 
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These are leading components of the fields (with the fall-off of order ry~^) 
with respect to the interpretation frame parallelly transported along a null geodesic 
approaching T, which has an affine parameter rj. The spherical angles 9, 4> describe 
the spatial direction on I of the null geodesic with respect to the reference tetrad 
(|29|l. q — cos^Qo + sin0 (cos^Tq + sin^So), and the parameter 9s is given by |3T|) . 
The coefficient is the part of ^'2 given by (fT4| that is proportional to ry"^, while 
the coefficient is the part of $1 given by (|18|) that is proportional to 77^^. Since 
= erj~^ near I [e — +1 for outgoing null geodesies while e = —1 for ingoing null 
geodesies), see equation (2.f9) in [fO], we immediately conclude that 



(1 — lUjpqY 



{m + in) + {e^ + g^) 



q ~ ap 
1 + i ujpq 



(34) 



^ e(l - ^Vg^) - 2gLupq + i [g(f - u;Vq^) + 2eu;pq] 
^* 2(1 + Lu^p^q^y ■ ^ ' 

Now, it only remains to evaluate the expressions (pij) . and ([55)1 on X. Using 
(|20p . that is by substituting g = — ap, all these functions become dependent only 
on p which specifies various points of scri for such (axially symmetric) spacetimes. 
Straightforward calculations lead to the following gravitational radiation component 

m\^^^gip)A{9,cl,,9s), (36) 



where 



with 



r( ^ 3 ;;^ a^P(p) + f(f + aW) 

gip) - 2 VW) A (1 + ,2^2^4)5/2 ' (^7) 

1 + a'^Lu^p'^ ' 1 + a'^uj'-^p'^ 



Vip) = + - 4«(e^ + ,^)^^^±^p + 4a^i^-t^p^ , (38) 



and ^(0, 0, 0s) is given by p3|) in which, using (|23|) 



Similarly, for the radiative component of the electromagnetic field we obtain 

= \s{p)Ad,'i^,Os), (40) 

where 

^(1 + a^w^p*)^ 

The function 7^(p) has the form (fTT|) with (l3])-([5]). These expressions give an explicit 
formula for the radiative components of the gravitational and electromagnetic fields 
in the complete family of Plebahski-Demiahski black hole spacetimes, characterized 
by the physical parameters to, a, a, I, e, g, A, in the "de Sitter-like" case A > 0. The 
amplitude of radiation Q [p) will be analyzed and discussed in sections 14.21 and 14.31 
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3. 3. The structure of radiation when A < and Qj > 

In the case when the cosmological constant A is negative and the acceleration 
parameter a is small, such that the condition (P5|) holds everywhere, all regions near 
the anti-de Sitter-like (timelike) conformal infinity are stationary and Qj > 0. 

To obtain the corresponding radiation components in these Plebahski-Demiahski 
spacetimes, we will first employ the special null frame k, 1, m, m given by (fT2|) for 
which the function Q is positive. The orthonormal tetrad tg, qs, fs, Ss associated with 
such principal null directions, introduced by (|27p. is thus 

ts = 7= [Or -Ujq da) , 

Q VQ 
Q. /— 
qs = —vQdq , 
Q 

n r- 

Ts = wVdp , 

Q 

Ss = — -j= ( ujp^dr +da). 

In the present case, it is convenient to choose the reference orthonormal tetrad 

to = --^{dr - ujq'^da) , 
9 VQ 



(42) 



(43) 



So = — {ujp'^dr + da) 



Obviously, tg = to, Sg = Sq. Using ((2T|) . qo = — n on X, i.e. qo is an outer normal 
to scri (see Fig. 1 and equation (4) in [11]), which corresponds to the choice cq = +1. 
Since —gap k^n^ = g^p l^n^ > in this case, the principal null direction k is outgoing 
with respect to X whereas the principal null direction 1 is ingoing, that is ei — ~ti, as 
in relation (29) in [11]. Moreover, 

qs = cosh" Vs qo + tanh V's fo , fs = - tanh t/Js qo + cosh" Vs i"o , (44) 

where 



sinhV's = -y'"^'^' cosh V's = y 1 + (45) 

which are relations of the form (31) and (10) therein. Thus, the magnitude of the 
radiative components of the gravitational and electromagnetic fields are given by 

l*4l = A^^i^^i^^''^'^^)' = ^^-A^{^,4>,i,:), (46) 

|?y| cosh V's n cosh V's 

where 

A\ {ip, V's) = (sinh ip + e sinh tps cos (f>)'^ + sinh^ t^s cosh^ V sin^0. (47) 

In these expressions, the pseudospherical parameters ip, (j) characterize the direction 
on timelike I of a given null geodesic with respect to the reference tetrad (1^5]) , namely 
t = cosh -0 to + sinh?/; (cos^Tq + sin(/)So). The parameter e indicates its orientation: 



Radiation generated by accelerating and rotating charged black holes 



11 



for e = +1 the direction is outgoing (pointing outside the spacetime) whereas for 
e = —1 it is an ingoing direction. 

The coefficient ^I^l^ in (I46p is again given by ([M]) . while is now the complex 
conjugate of (|35p . Therefore, the radiative components of the gravitational and 
electromagnetic fields in this case can be written as 

\%\ ^^\g{p)\Ai{^A,A), (48) 

\v\ 

^ ^\S{p)\Ai{i;,ct>,i^s), (49) 

where the functions G{p) and £{p) have been introduced in (|37p and (|¥T|l . The 
directional structure Aiiip, (j),ijjs) is given by expression (|47p in which 

sinh^ - > (50) 

-f(H-a2wV)-a27'(p) ^ ' 

determines the orientation of the principal null directions with respect to the conformal 
infinity I. These are explicit formulae for the radiative components of the fields in 
the family of Plebahski-Demiahski black hole spacetimes in the "anti-de Sitter-like" 
case A < when the acceleration parameter a is small. In such a situation, Q > 
everywhere near X, on which the principal null direction k is outgoing while the 
principal null direction 1 is ingoing. In the subcase uj = 0, the results obtained agree 
with those derived in [7], i.e. for the C-metric with < which represents a single 
non-rotating but accelerating black hole in anti-de Sitter universe. 

3.4- The structure of radiation when A < and Qx < 

When the cosmological constant A is negative but the acceleration parameter a is 
large, there are also regions near the anti-de Sitter-like (timelike) conformal infinity in 
which the condition holds. In such regions where Qj < 0, the spacetime is not 
stationary and the asymptotic structure of radiation is different from that described 
in previous section \3l3[ 

To determine the corresponding radiation components we will employ the 
algebraically special null frame closely related to p^ . namely l,k, — rn, — m. The 
associated orthonormal tetrad tg, qs, Tg, Sg is thus 

n , 

ts = -\/-Qa,, 

{dr - ujq^d„) , 

(51) 
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cf. relations (|771) and A convenient reference orthonormal tetrad is 



n -VQ 



^VQ + a^V (52) 

:{dr-UJq^d^), 



So = j= ( Ujp dr + Oa) ■ 

As in the previous section, using (|2ip and we obtain Qo = — n on X , so that 
Co = +1. However, Q < which means that qo is now the inner normal to scri. Also, 
Qap k"n^ = gap l^n^ < in this case, so that both principal null directions k and 1 
are ingoing with respect to 2 (and ei = £2, as in relation (40) in [11]). Moreover, 
Qs = To, Sg = So, and 

ts — cothV^s to + sinh^"'^0s ^o, — ~ sinlC^iJjs to — coth-^s Qo, (53) 

where 



cosh-fAs = •y'"^!^ ' sinhV's = ya^^^ - 1 ■ (54) 

These are particular exemplifications of relations (42), (10) of [11]. Consequently, the 
radiative components of the fields are given by expressions 

m\^T^^^M^A,A), m^^^^^^A2{^A,A), (55) 

\i]\ smh Vs V smh V's 

where 

A2{ip,4>,i^s) = (coshes + e cosh f/;)^ + sinh^ sinh^T/isin^^i . (56) 

Again, the pseudospherical parameters tp, (j) characterize the direction on timelike X 
of a given null geodesic with respect to the reference tetrad (|52|) . e = ±1 denotes its 
orientation, the coeflrcient is given by p4|) . while is given by ([35| multiplied 
by the overall factor (—1). 

The radiative components of gravitational and electromagnetic fields in this case 
are thus 

\%\ ^^\g{p)\A2{^A,A), (57) 

\v\ 

= 4l^(p)|A(V',0,V's), (58) 

where the functions g{p) and £{p), now with an opposite sign, have been defined in 
(|37p. (|4T|) . The directional structure A2{'4',(f>,tps) is now given by expression ((56| in 
which, using ^ and (|23ll . 

a^r{p) 



'^o^^ ^« ^2p(p) + A(i + a^Lo^pi) ' 
a^^(p) + f(l + a^^V) 



These are explicit formulae for the radiative components of the fields in the family 
of Plebanski-Demiahski black hole spacetimes in the "anti-de Sitter-like" case A < 
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when the acceleration parameter a is large and Q < 0. More precisely, for large a there 
are both stationary (Q > 0) and non-stationary (Q < 0) regions near the conformal 
infinity T. In the present case Q < 0, so that both the principal null directions k 
and 1 are ingoing on X. Note that for a; = 0, the above results reduce to 

those derived in [7], i.e. for the C-metric with > which describes a pair of 
non-rotating, uniformly accelerating black holes in anti-de Sitter universe. 

4. Discussion of radiative properties of the Plebanski Demianski black 
hole spacetimes 

In the preceding section [3] we have derived explicit forms of radiative components of 
gravitational and electromagnetic fields for a completely general family of Plebanski- 
Demiahski black hole spacetimes with a cosmological constant and a twist. The 
asymptotic structure of such fields has been described by their specific components 
evaluated in the interpretation frame that is parallelly transported along null geodesies 
approaching T. We have denoted them as 1^*4 1 and |$2L respectively. 

When the cosmological constant is positive (A > 0) the conformal infinity T has 
a de Sitter-like character (i.e., its normal is everywhere timelike) and the asymptotic 
components of the fields are given by expressions ([36| , (j40| . In the anti-de Sitter- like 
case (A < 0) the normal to scri is spacelike, and it is necessary to distinguish two 
cases. For a small acceleration a of the black holes, more precisely in the points of T 
characterized by the coordinate p where a^{'P{p) -\- ^w^p*) < — , the asymptotic 
fields are given by (pS)) , (|49p , while for a large acceleration there are also regions on I 
where a^{V(p) + f w^p^) > , in which the fields are given by (H?]), 

All these radiative components asymptotically decay as where rj is an affine 
parameter along the null geodesic (and jryj ^ oo as the geodesic approaches I). 

4.I. The directional structure of radiation on X 

When A > 0, the directional structure of radiation at any point at conformal infinity 
is fully determined by the function A{d, cj), 9s) given in (f33|) . In this case, the spherical 
angles 0, (f> describe the spatial direction q = cos qo + sin 9 (cos 4>ro + sin (pSo) on T 
of a given null geodesic with respect to the reference orthonormal tetrad (p9)) . The 
parameter 6s is explicitly determined by ((39)) — let us recall that the variable p 
therein determines the specific point on I where the null geodesic "terminates" . Such 
a directional structure of radiation (which is visualized on Fig. 3 of [11]), is typical for 
all type D spacetimes that admit a spacelike conformal infinity. 

When A < and acceleration a is small, the directional structure of 
radiation on a timelike conformal infinity is described by the function Aiiip, (j),ips) 
introduced in ([47|) . The pseudospherical parameters ip, (j) characterize the 
direction on X of a given null geodesic with respect to the reference tetrad ([43]) . 
t = cosh-0to -I- sinh?/; (cos0ro -I- sin(/)So), and Vs given by (l50|) determines the 
orientation of the principal null directions with respect to X (here k is outgoing while 1 
is ingoing since Q > everywhere near X). However, when the acceleration a is large, 
there are also non-stationary regions near timelike conformal infinity where Q < 
(and both k and 1 are outgoing), in which case the asymptotic directional structure 
of radiation at a given point on X is described by A2{'>p^ 4>,^s)- This is given by (|56|) . 
with the parameter ^s determined by (|59p . Again, these directional structures (shown 
on Fig. 5 of [11]) are characteristic for type D spacetimes with a timelike X. 



Radiation generated by accelerating and rotating charged black holes 



14 



4-2. The amplitude of radiation on 2 

The main result of the present work is the derivation of exphcit asymptotic radiative 
component of the gravitational field in the complete family of Plebahski-Demiahski 
black hole spacetimes with a possible twist and a cosmological constant A. This is 
uniquely determined by the amplitude G{p) given by ([37|) . The analogous radiative 
component of the electromagnetic field is given by the amplitude £{p), see (|4T|). 
Interestingly, the functions \G\ and \£\ are the same for all three possibilities discussed 
above, namely A > and A < with either Q > or Q < 0. 

In the remaining part of this contribution we will thus concentrate on description 
and analysis of such amplitudes. In particular, we will investigate and visualize these 
functions which uniquely express the dependence of the magnitude of gravitational 
and electromagnetic radiation on the parameters that characterize various physical 
properties of the black holes, namely their mass m, electric and magnetic charges 
e and g, NUT parameter rotational parameter a, and acceleration a. Within 
the framework of rigorous Einstein theory, we may thus study specific properties of 
radiation generated by localized sources of the black hole type, as described by the 
metric ([T0| . ([TT1) . 

The amplitude of gravitational radiation Q (p) is given by expression (I37p in which 
the specific function takes the form ([55)1 . The parameter p, together with the 

angular coordinate a — -~{uj/a) (j), where € [0, 2tt), labels all possible points on the 
(axially symmetric) conformal infinity X. As explained in section [21 the coordinate 
p must be restricted to a particular range between appropriate roots V = which 
correspond to poles of the black holes. This is automatically achieved by assuming 
the parameterization p — (l/uj) + {a/oj)cos0 with 9 S [0, tt], see transformation 
Simultaneously, it is necessary to express the (unphysical) coefficients e, n and k in 
(fTTjl using the expressions ©-(IS]). 

The resulting general expression for G{p{d)) is involved. For simplicity, we will 
thus restrict our attention only to the case, in which the NUT parameter I vanishes. 
Setting I = 0, the black hole spacetimes are free from torsion singularities, regions with 
closed timelike curves are absent, and the spacetimes are physically more realistic. 
With this assumption, the twist parameter ld can be chosen to be equal to a, so that 

p^msO. (60) 

The complete range 9 £ [0, tt] thus corresponds to p E [—1, +1]. The values p = ±1 
represent poles of the black hole spacetimes, while the special value p — describes 
the equatorial plane 9 — ^. 

In such a case, expressions ([S])-® simplify considerably to 

fc = l, n = -aam, e = 1 - a^(e^ + g^) - (a^ + |)a^, (61) 

the function V{p), given by (fTTjl . becomes 

V{p) = (1 - p^)(l - 2amp+ [a^{e^ + g^) + (a^ + ^)a^ ] , (62) 

and the function T>{p), introduced in ([55]) . takes the form 



V{p) = m\l + a^a') - Aam{e^ + g^)\^^^p + Aa^ f/^f^'y . (63) 

l + a^a^p^ 1 + a-^a^p* 

The amplitude ((37)) is therefore 



a \ 3 a^7^(p) + f(l + AV) 
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The above expressions explicitly describe the amplitude of gravitational 

radiation in the spacetimes with a cosmological constant A which represent uniformly 
accelerating and rotating charged black holes, that are characterized by the parameters 
a, a, e, g and m, respectively. In particular, in the equatorial plane given by p = 0, 
the amplitude simplifies considerably to 

g{0) ^^rny/l + a^a^ ( 1 + fa^) . (65) 

In the case of uncharged black holes (e = = g), the amplitude of gravitation 
((64)) reduces to 

g{p) - 2 m Vl + a^a^ a.. ^ \ ^ 4.5/2 , (66) 

■3 (1 + Q!^a^p*)°/^ 

where 

V{p) = (1 - p^) (l - 2amp + (a^ + A)a2 ^g^^ 

These expressions describe the amplitude of gravitational radiation generated by Kerr- 
like black holes (of mass m and rotation parameter a) which uniformly accelerate (as 
determined by the acceleration parameter a) in the de Sitter or anti-de Sitter universe 
with a cosmological constant A. 

For a vanishing rotation of the black holes, a = 0, this further simplifies to 

Gip) ^lm(l + j^a\l~p')il-2amp)^ , (68) 

which corresponds to the uncharged C-metric, while for a vanishing acceleration, 
a = 0, we immediately obtain 

0s =0 when A > 0, (69) 
tps^O when A < 0, 

which corresponds to the Kerr-(anti-)de Sitter black holes. In this case, the amplitude 
of gravitational radiation Q is constant (proportional to the mass parameter m of 
the black hole), independent of the rotation parameter a and the specific value of 
the cosmological constant A. Moreover, it is also independent of the position p on 
scri T. The directional structure of radiation on X is given by A(B^ = sin^ B or 
^i('0) = sinh^?/), if A is positive or negative, respectively. 



Jf..3. Visualisation of the results for A > 

In this case, for a generic choice of the parameters of rotating charged black holes 
which accelerate in de Sitter universe, the function G{p) determined by the expression 
(|64p is plotted in figure[TJ As a typical example, the values are taken as A = 1, m = 1, 
a = 0.5, e = 0.5 = g, and a — 0, 0.1, 0.2, 0.3, 0.4, 0.5. It can be seen that the amplitude 
of gravitation radiation has a single maximum around the value p ~ —0.5 (i.e. in the 
"southern hemisphere" of I ) and decreases monotonously to a minimal value at the 
north pole p — I- The greater the acceleration a, the bigger the generated radiation. 

Similar results are obtained for uncharged rotating black holes (e = = 5) , see 
figure [2] plotting and non-rotating black holes {a — 0), see figure [3] for the case 
of an uncharged C-metric. While the dependence on the acceleration a is quite 
strong, the dependence on the charge parameters e, g is weaker, see figure SI and 
the dependence on the rotation parameter a is much weaker, see figure [5j 
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Figure 1. The amplitude of gravitational radiation G{p) on a spacelike I, in 
which case Qx < for any point p of X. As a typical example, it is assumed that 
A = 1, m = l, a = 0.5, e = 0.5 = g, while the acceleration of black holes is taken 
in the interval a G [0,0.5]. The horizontal line at Q = 1.5 represents the uniform 
amplitude for non-accelerated Kerr-Newman-de-Sitter spacetime with a = 0. 





Figure 3. The amplitude G{p) on a spacelike X for A = 1, m = 1, a = 0, 
e = = g, and a £ [0,0.49]. This represents radiation of non-rotating uncharged 
black holes which accelerate in de-Sitter spacetime, described by the C-metric. 
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Figure 4. The amplitude G{p) on a spacelike X for A = 1, m = 1, a = 0.5, 
= 0.3, and e = g£ [0,0.5]. This shows the dependence of gravitational 
radiation, generated by accelerating rotating black holes, on their charges. 
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Figure 5. The amplitude G{p) on a spacelike X for A = 1, m = 1, e = 0.5 = g, 
a = 0.3, and a S [0, 0.5]. The dependence of gravitational radiation, generated by 
accelerating charged black holes, on their rotation a is very weak. 

4-4- Visualisation of the results for A < and Qj > 

When the cosmological constant A is negative and the acceleration parameter a is 
small enough, such that the condition (P5)l holds everywhere on I, all points p on 
the anti-de Sitter-like (timelikc) conformal infinity are in the stationary region since 
Qx > 0. For a typical choice of the parameters of accelerating rotating charged black 
holes the function Gip) is plotted in figure [HI As an example, the values are taken 
as A = -1, TO = 1, a = 0.5, e = 0.5 = g, and a = 0, 0.1, 0.2, 0.3, 0.4, 0.5. It can be 
seen that for very small values of a, the amplitude of gravitation radiation decreases 
uniformly from the south pole {p — —1) to the north pole {p = +1), while for larger 
accelerations there is a significant local minimum in the "southern hemisphere" . 

Analogous results are obtained for uncharged rotating black holes {e = = g), see 
figure [71 and non-rotating black holes (a = 0), see figure [8] which corresponds to the 
case of an uncharged C-metric. As in the previous case A > 0, the dependence on the 
acceleration a is generally strong, the dependence on the charges e, g is weaker, see 
figure [9l and the dependence on the rotation a is considerably weaker, see figure fTOl 
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Figure 6. The amplitude of gravitational radiation G{p) on a timelike X in the 
case when Qx > 0. As a typical example, it is assumed that A = —1, m = 1, 
a = 0.5, e = 0.5 = g, while the acceleration of black holes is taken in the interval 
o S [0,0.5]. The horizontal line a,t Q = 1.5 represents the uniform amplitude for 
non-accelerated Kerr— Newman-anti-de-Sitter spacetime with a = 0. 
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Figure 7. The amplitude Q{p) on a timelike I and Qx > 0, for A = —1, m = 1, 
a = 0.5, e = = g, and a G [0, 0.48]. This describes gravitational radiation of 
uncharged rotating black hole which accelerate in anti-de-Sitter spacetime. 
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Figure 8. The amplitude on a timelike X and Qx > 0, for A = —1, m = 1, 
a = 0, e = = g, a S [0, 0.49]. This describes radiation of non-rotating uncharged 
black holes which accelerate in anti-de-Sitter spacetime (given by the C-metric). 



Radiation generated by accelerating and rotating charged black holes 



19 



a=e=g=0.0 

2.0 - e=g=0.1 

e=g=0.2 




Figure 9. The amplitude G{p) on a timelike X and Qx > 0, for A = —1, m = 1, 
a = 0.5, a = 0.3, and e = g G [0, 0.5] which shows the dependence of gravitational 
radiation, generated by accelerating rotating black holes, on their charges. 

2.2 
2.0 
1.8 
"1.6 



1.4 
1.2 



1.0 

-1 -0.5 0.5 1 

P 

Figure 10. The amplitude Q{p) on a timelike X and Qj > 0, for A= — 1, 
m = 1, e = 0.5 = g, a = 0.3, and a £ [0, 0.5]. This shows that the dependence 
of gravitational radiation on rotation of these black holes is very weak. 



4-5. Visualisation of the results for A < and Qj < 

When the cosmological constant A is negative but the acceleration parameter a is 
large, the condition (|25p for Qj > is not satisfied everywhere on the anti-de Sitter- 
like (timelike) conformal infinity X. In the regions (i.e. for appropriate ranges of p) 
where Qj > 0, the amplitude of gravitational radiation G{p) still behaves qualitatively 
in the same way as described in section H^ However, now there are also non-stationary 
regions of I in which the condition ((26l) holds so that Qj < 0. (For more details about 
the structure of the anti-dc Sitter I for large values of a see [7,20].) 

For a typical such case A = —0.7, m — 1, a = 0.2, e ~ 0.5 = g, and a ~ 0.45, 
0.48,0.51,0.54,0.57, the function \Gip)\ is plotted in figure fTTl Admissible ranges of 
the coordinate p and of the acceleration parameter a result from the necessity to satisfy 
both the conditions V > and Qx < 0. There is a single maximum localized on the 
southern hemisphere, approximately in the centre of the corresponding interval of p. 
It can be seen from figures fT2l and [T3l that for uncharged {e — — g) and non-rotating 
(a = 0) black holes the behaviour of the function \G{p)\ is basically the same. 
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Figure 11. The amplitude of gravitational radiation on a timelike I in the 

regions where Qj < 0. It is assumed that A = —0.7, m = 1, a = 0.2, e = 0.5 = g, 
while the acceleration of the black holes is taken in the interval a G [0.45,0.57]. 
For completeness, we also plot in the stationary regions where Qx > by 

the corresponding grey lines. 




Figure 12. The amplitude on a part of timelike X where Qx < 0, for 

A = —0.7, m = 1, a = 0.2, e = = g, and a £ [0.45, 0.50]. It describes radiation 
of uncharged rotating black holes accelerating in anti-de-Sitter spacetime. 




Figure 13. The amplitude |5(p)| on a part of timolike X where Qx < 0, for 
A = —0.7, m = l, a = 0, e = = g, a£ [0.45, 0.49], describing radiation of non- 
rotating uncharged black holes accelerating in anti-de-Sitter space (the C-metric). 
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Again, generically the dependence on the acceleration a is strong, the dependence 
on the charges e, g is weaker (figure fn)) . and the dependence on the rotation a is almost 
negligible (figure [T5|) . 




Figure 14. The amplitude of gravitational radiation |(7(p)| on a part of timelike 
I where Qj < 0, for A = -0.7, m = 1, a = 0.2, a = 0.5, and e = g e [0, 0.5]. This 
shows the dependence on charges of the accelerating black holes. 
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Figure 15. The amplitude of gravitational radiation |(7(p)| on a part of timelike 
I where Qx < 0, for A = -0.7, m = 1, e = 0.5 = g, a = 0.5, and a e [0, 0.5]. This 
shows that the dependence on rotation of the accelerating black holes is very weak. 




4-6. Electromagnetic radiation 

The amplitude of electromagnetic radiation £{p) is determined by expression ((4T|) . in 
which the coordinate p labels all possible points on the (axially symmetric) conformal 
infinity T. For a generic choice of specific parameters of accelerating and rotating 
black holes, this function is visualised as in figures [TBHT51 Figure [11] applies to the 
case when A > while the remaining two figures apply when A < 0. Specifically, 
figure [T7| describes the function £{p) in those regions of X where Qx > 0, and figure [TSl 
in the points where Qx < 0. It can be seen that the behaviour of the electromagnetic 
amplitude £{p) is qualitatively similar to the behaviour of the gravitational amplitude 
G {p) , cf . figures m [5] and [111 respectively. 




Figure 16. The amplitude of electromagnetic radiation £(p) on a spacelike X, 
in which case Qx < for any point p of I. Here it is assumed that A = 1, 
m = 1, a = 0.5, e = 0.5 = g, while the acceleration of the black holes is taken in 
the interval a £ [0,0.5]. The horizontal line at £ = 0.125 represents the uniform 
amplitude for non-accelerated Kerr-Newman-de-Sitter spacetime with a = 0. 
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Figure 17. The amplitude of electromagnetic radiation £(p) on a timelike X in 
the points where Qx > 0. Here A = —1, m = 1, a = 0.5, e = 0.5 = g, while the 
acceleration of the black holes is taken in the interval a € [0,0.5]. 




Figure 18. The amplitude of electromagnetic radiation £{p) on a timelike X in 
the points where < 0. Here A = —0.7, m = 1, a = 0.2, e = 0.5 = g, while the 
acceleration of the black holes is taken in the interval a £ [0.45,0.57]. 
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It should also be mentioned that the amplitude and directional structure of 
radiation, as described by expression (HD|) . gives (47r-times) the magnitude of the 
Poynting vector in the parallelly transported interpretation frame. Moreover, it 
exactly reduces to results previously obtained for non-rotating black holes [6, 7], 
and further for test fields generated by uniformly accelerated charges in de Sitter 
spacetime [4, 5] in the case when the mass of the sources becomes negligible. 

5. Conclusions 

In this contribution we analyzed in detail the asymptotic behaviour of gravitational 
and electromagnetic fields of black hole spacetimes which belong to the Plebahski- 
Demiahski family of solutions. The amplitude and directional structure of radiation 
on conformal infinity X were evaluated in cases when the cosmological constant is non- 
vanishing, so that I have either de Sitter-like or anti-de Sitter-like character. Although 
the accelerated systems studied do not seem to be astrophysically realistic (because 
their acceleration is caused by cosmic strings or struts attached to them along the axis 
of symmetry), the results obtained elucidate some principal geometric and physical 
properties of this class of exact radiative spacetimes. 

In particular, we found explicit relations between the parameters that characterize 
the sources, namely their mass, electric and magnetic charges, NUT parameter, 
rotational parameter, and acceleration, and the asymptotic radiative properties. For 
example, it was demonstrated that radiation, described by specific components of 
the Weyl and Maxwell tensors in the parallelly transported interpretation frame near 
conformal infinity, is not distributed uniformly on T when the black holes accelerate. 
Also, both the gravitational amplitude G{p) and the electromagnetic amplitude £{p) of 
radiation on I depend strongly on the acceleration of the black holes, weakly on their 
charges, and very weakly on their rotation. Since previously only non-rotating black 
holes (described by the C-metric) were studied, this is a new and rather surprising 
result. 

As already mentioned above, there are fundamental differences concerning the 
asymptotic structure of radiation in spacetimes with and without a cosmological 
constant. Since the geometric nature of the corresponding conformal infinities is 
different it is not straightforward to perform the limit A — > 0. Indeed, the amplitudes 
G{p) and E{p), given by ([77)1 and (|4ip . diverge in this limit. Nevertheless, in the case 
A = it is possible to evaluate directly the dominant (radiative) components of the 
fields in an appropriate interpretation frame. It follows from (1251) that Q < near the 
conformal infinity (which now has a null character), so that the region must be non- 
stationary. Using then the reference frame to = ts,qo — rs,ro ~ qs,So = Sg, where 
ts,qs,rs,Ss are given by ([28|) . it turns out that the gravitational and electromagnetic 
fields behave asymptotically as l^i^l « 3|*|J \t]\~^ and ~ V2|<i>fJ |??r\ where 
the coefRcients are given by (p4|) and (|35|) . i.e., |'I'|^| — ^J'D{p) (1 + a^u'^p^)^^ and 
|*!*| = \\/^Tf{l a2wV)"\ with V(p) defined by Compared to ^ and 
(|40p. there is obviously no directional structure when A = (the fields are the same 
for all null geodesies approaching a given point at null conformal infinity). Also, with 
the above normalisation, the dependence on p (which labels all possible points at X) is 
simpler and the amplitudes of gravitational and electromagnetic radiation are finite. 
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